Electronic structure and magnetic anisotropy for nickel-based molecular magnets 
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Recent magnetic measurements on tetra-nickel molecular magnets [Ni(hmp)(ROH)Cl]4, where 
R=CH3, CH2CH3, or (CH2)2C(CH3)3 and hmp~ is the monoanion of 2-hydroxymethylpyridine, 
revealed a strong exchange bias prior to the external magnetic field reversal as well as anomalies 
in electron paramagnetic resonance peaks at low temperatures. To understand the exchange bias 
and observed anomalies, we calculate the electronic structure and magnetic properties for the N14 
molecules with the three different ligands, employing density-functional theory. Considering the 
optimized structure with possible collinear spin configurations, we determine a total spin of the 
lowest-energy state to be S = 0, which does not agree with experiment. We also calculate magnetic 
anisotropy barriers for all three types of M4 molecules to be in the range of 4-6 K. 

PACS numbers: 75.50.Xx, 73.22.-f, 75.30.Gw, 71.15.Mb 
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Molecular magnets have received much attention in 
the past decade because of observed quantum features 
at the mesoscopic level and potential utilization as 
ultra-high density information storage devices or mag- 
netic sensors. The first prototype molecular magnet is 
Mni 2 Oi2(CH3COO)i 6 (H 2 0)4 (hereafter Mni 2 -acetate),± 
which has a total ground-state spin of S = 10 and mag- 
netic anisotropy barrier of 65 K^ 3 Magnetic hystere- 
sis loop measurements on Mni2-acetate showed quan- 
tum tunneling of magnetic moment through the magnetic 
anisotropy barrier when the external magnetic field was 
applied along the easy axis^ 

Recently molecular magnets [Ni(hmp)(ROH)Cl]4 
(hereafter N14) were synthesized^ where R=CH 3 , 
CH 2 CH 3 , or (CH 2 ) 2 C(CH 3 ) 3 (abbreviated to Me, Et, 
and tBuEt, respectively). They revealed qualitatively 
different tunneling features from the typical molecular 
magnets such as Mni 2 -acetate. For N14, the tunneling 
was observed prior to magnetic field reversal but no 
tunneling in zero magnetic field^ The similar behavior 
has been observed in a dimeric form of the molecular 
magnet [M^OgCU^CEt^NCsHg^^ (hereafter 
Mn.4)m& Additionally, electron paramagnetic resonance 
(EPR) spectra on N14 exhibited unusual double sets 
of low-temperature peaks corresponding to somewhat 
different values of magnetic anisotropy barriers^ This 
anomaly in the EPR spectra has never been reported 
from any other kinds of molecular magnets including 
the Mn4 dimer. Each N14 molecule comprises four 
Ni 2+ (Si = 1) ions located at the corners of a slightly 
distorted cube and coupled through oxygen anions with 
S4 symmetry shown in Fig. ^ Unlike Mn 12 -acetate 
and Mn4 dimer, the N14 molecular magnets are not 
well separated from each other so that there are four 
interacting nearest neighboring molecules ^ The inter- 
molecular interactions can be tailored by different sizes 
of ligands such as MeOH, EtOH, and tBuEtOH. Among 
the three different types of Ni 4 , the MeOH and EtOH 



complexes consist of two sublattices, each of which is an 
inter-penetrating diamond latticed 




FIG. 1: Geometry of molecule magnet [Ni(hmp)(MeOH)Cl] 4 . 
For simplicity, hydrogen atoms are not shown. 

In this paper, as an attempt to better understand the 
observed tunneling aspects and EPR spectra, preliminary 
first-principle calculations on the Ni 4 molecular magnets 
will be presented with the three different ligands using 
density- functional theory (DFT). Geometries of the Ni 4 
molecular magnets with the different ligands will be re- 
laxed and electronic density of states will be discussed 
along with ground-state magnetic moments and magnetic 
anisotropy. All of our calculated results will be compared 
with available experimental data. 

Density-functional calculations 8 were performed with 
spin-polarized all-electron Gaussian-orbital basis sets and 
the generalized-gradient approximation (GGA)i^For this 
purpose, the Naval Research Laboratory Molecular Or- 
bital Library (NRLMOL)A2i is used. Compared to Mn i2 - 
acetate and M114 dimers, the Ni 4 molecular magnets 
need much finer mesh for faster convergence of total en- 
ergy. Five different Ni 4 molecules are taken into account 
in density-functional calculations: a tBuEtOH complex 
and two crystallographically inequivalent molecules of 
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TABLE I: Calculated energies (in units of eV) of different 
collinear spin configurations relative to the lowest-energy con- 
figuration (M s = 0) for the three types of N14. Here MeOH.01 
represents a molecule belonging to the first sublattice of the 
MeOH complex. 





E{M S = 4) 


E(M S = 2) 


E(A/ S = 0) 


MeOH.01 


0.0308 


0.0023 


0.0000 


EtOH.Of 


0.0289 


0.0042 


0.0000 


tBuEtOH 


0.0283 


0.0056 


0.0000 



TABLE II: Magnetic moments captured by spheres around 
ions in units of fis for M s = 4 of the EtOH complex and 
for M s = 10 of the Mni2-acetate. For N14, O(b) denotes a 
bridging oxygen and 0(1) represents a dangling oxygen from 
the cube. For Mm 2 -acetate, 0(1) [0(2)] bridges Mn 4+ (1) 
[Mn 3+ (2)] ions. 



Ni 


CI 


O(b) 


0(1) 


N 




1.57 


0.053 


0.166 


0.043 


0.049 




Mn 4 +(1) 


Mn 3+ (2) 


Mn 3+ (3) 


0(1) 


0(2) 


0(3) 


-2.58 


3.62 


3.55 


0.012 


-0.024 


-0.029 



the MeOH and EtOH complexes. For each molecule, 
a relaxed geometry is found starting from a ferromag- 
netic M s — 4 spin configuration until forces exerted on 
any atom become relatively small (where 2M S is mag- 
netic moment). Although the moment is not fixed during 
self-consistent iterations, it converges to M s = 4. With 
the relaxed geometry two other collinear spin configu- 
rations, M s = 2 and M s = 0, are considered as start- 
ing magnetic moments to calculate energies and com- 
pare to the energy of M s =4. As shown in Table [I] 
for the MeOH, EtOH, and tBuEtOH complexes, the 
M s = configuration has the lowest energy. This is 
in contrast to the experimental result*^ which exhib- 
ited a ferromagnetic coupling S = 4 as the ground state 
for all five Ni 4 molecules. Even though not shown in 
Table Q] there are no differences in total moments of 
the ground states for the two inequivalent molecules of 
the MeOH and EtOH complexes. An experiment on a 
tetra-nickel-based molecular magnet within a M012 cage, 
[Moi 2 03o(^2-OH)ioH 2 {Ni(H 2 0)2}4], however, revealed 
that the ground state has a total spin of S = 0>ii The 
experimental finding of the ground-state spin was con- 
firmed by density-functional calculations carried out by 
Postnikov et alM> In this system, two Ni ions are bridged 
via O-M02-O, while in our system Ni ions are bridged via 
oxygen anions only. 

To investigate what may cause the discrepancy be- 
tween theory and experiment, we calculate the magnetic 
moment captured by each atom, the electronic density of 
states projected on individual atom, and the spin-flip en- 
ergy gaps for the M s = 4 configuration. Within a sphere 
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FIG. 2: Electronic majority- and minority-spin density of 
states (DOS) projected onto Ni(3d), Ni(4s), Cl(2p), 0(2p), 
and N(2p) of [Ni(hmp)(EtOH)Cl] 4 . All projected DOS have 
the same scale except for Ni(3d) DOS. The vertical line de- 
notes the Fermi level. 



of 2.3 bohr radius around a Ni ion, a moment of 1.6/is 
is found. The second dominant contribution to the total 
moment arises from oxygen anions bridging Ni ions in 
the cube. A moment of 0.17 /is is found within a sphere 
of 1.2 bohr radius around one of the oxygen anions. This 
is approximately 10% of the captured moment around 
the Ni ion. Compared to a ferrimagnetic structure of 
the Mni2-acetate and M114 monomer ^iiiiiiiSi somewhat 
larger moments are found around the oxygen anions in 
Ni 4 (Table [HJ . Interestingly, for the M s = configu- 
ration, a bridging O anion in M4 bears a moment of 
0.073 fiB within a sphere of 1.2 bohr radius. From the 
calculated electronic density of states (Fig. |5Jl projected 
onto Ni 3d orbitals for M s = 4, we find that most of the 
electron spin density is localized at the Ni sites. Notice 
that some of the oxygen 2p orbitals for minority spin are 
not occupied, although all majority-spin O 2p orbitals 
are occupied. This might imply a possible small leakage 
of spin density to the O sites. However, more scruti- 
nized studies are required to conclude that this, in fact, 
attributes to the antifcrromagnetic coupling in the cal- 
culated lowest-energy state. Another intriguing point is 
that as shown in TableQJ collinear spin excitations for Ni 4 
have an order of magnitude lower energies than those for 
the Mni2-acetateM and Mn 4 monomeipi£*i&. A majority- 
minority spin-flip gap ranges from 1.4 eV to 1.6 eV and 
a minority-majority spin-flip gap varies between 1.8 eV 
and 2.1 eV depending on ligands as shown in Table ITTT1 It 
seems that the M s = 4 configuration appears to be mag- 
netically stable since the spin-flip gaps are much larger 
than thermal energy. 

An expectation value of the orbital angular momen- 
tum in the ground state of molecular magnets is typi- 
cally quenched so that the ground-state spin manifold 
has a (2S + l)-fold degeneracy. The degeneracy can be 
lifted by spin-orbit coupling, which leads to magnetic 
anisotropy in the system. A generalized second-order 



3 



TABLE III: Majority LUMO (lowest unoccupied molecular 
orbital)- HOMO (highest occupied molecular orbital) gaps, 
minority LUMO-majority HOMO gaps, majority LUMO- 
minority HOMO gaps, minority HOMO-LUMO gaps in units 
of eV, and magnetic anisotropy barriers (MAB) in units of 
kelvin for M a = 4 of the three types of Ni.4. 





Mj L-H 


Mn L-Mj H 


Mj L-Mn H 


Mn L-H 


MAB 


MeOH.Ol 


1.74 


1.38 


1.83 


1.47 


4.4 


EtOH.Ol 


2.02 


1.47 


2.18 


1.62 


5.9 


tBuEtOH 


1.97 


1.57 


2.12 


1.71 


5.4 



single-spin Hamiltonian is given by 

H = -DSl + E{S 2 X - S 2 y ) (1) 

where z is the easy axis, D is the uniaxial anisotropy 
parameter, and E is the transverse anisotropy parame- 
ter. To calculate the energy shift of the ground state, 
we consider the spin-orbit coupling as a small perturba- 
tion and use the computed single-electron orbitals and 
orbital energies^ for the relaxed M s = 4 geometry. Be- 
cause of S4 symmetry of the M4 molecule, the first term 
in Eq. only survives. Calculated magnetic anisotropy 
barriers (— DM 2 ) for all five different N14 molecules are 
in the range of 4-6 K so that D = 0.25 - 0.37 K. See Ta- 
ble ^Qj Measured EPR spectra on the five Ni 4 molecules 
revealed that values of D are in the range of 0.72-1.03 K 
and that the two inequivalent molecules of the MeOH 
and EtOH complexes have considerably different values 
of Dm- The calculated values of D are merely 36% of the 



experimentally extracted values. Numerical uncertainties 
in density-functional calculations prevent from differenti- 
ating between the anisotropy barriers of the two inequiv- 
alent molecules of the MeOH and EtOH complexes or 
between the barriers of the N14 molecules with different 
ligands. A single-ion anisotropy projectecU^ii onto a Ni 
ion in the MeOH complex provides local D = 1.52 K 
and local E = 0.47 K. The local easy axis of the Ni ion 
is found to be tilted by 10.5° from the global easy axis. 
Due to lack of Jahn- Teller distortion, a small value of the 
local D is theoretically expected for Ni 2+ ions. 

In summary, we have considered the N14 molecular 
magnets with three different ligands in volume and car- 
ried out first-principle calculations of the electronic and 
magnetic structure. We found that the calculated lowest- 
energy state has a total spin of S = 0, although ex- 
perimental data indicated S = 4. This is quite un- 
usual because density-functional theory, so far, pre- 
dicted/confirmed correct ground-state spins for a variety 
of molecular magnets such as the Mni^-acetat o 1 ?! 14 , M114 
monome»i£ii£, C04 molecule^, and Fe4 molecule^. Al- 
though the reason has not yet been clarified, our prelimi- 
nary studies tentatively suggested a plausible incomplete 
localization of spin density on Ni ions. There were no 
prominent effects of the ligands in different volume on 
the electronic structure and magnetic anisotropy of the 
N14 molecules. 
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